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We establish the physical process version of first law by studying small perturbations of a station- 
ary black hole with regular bifurcation surface in Einstein-Gauss-Bonnet (EGB) gravity. Our result 
shows that when the stationary black hole is perturbed by a matter stress energy tensor and finally 
settles down to a new stationary state, the Wald entropy increases as long as the matter satisfies 
null energy condition. 
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The striking similarity of the laws of black hole me- 
chanics with thermodynamics was first established in 
case of general relativity (GR) [![. It is interesting to 
explore whether this analogy is a peculiar property of 
GR or a robust feature of any generally covariant theory 
of gravity. 

The zero'th law, which ensures the constancy of surface 
gravity is valid for any stationary Killing horizon with a 
regular bifurcation surface Q, irrespective of the gravi- 
tational dynamics. 

The equilibrium state version of first law is established 
by Wald and collaborators [!, H| for any arbitrary diffco- 
morphisom invariant theory of gravity. Comparing with 
the first law of thermodynamics, the entropy of the black 
hole can be expressed as a local geometric quantity inte- 
grated over a space-like cross section of the horizon and 
is associated with the Noether charge of Killing isometry 
that generates the horizon. 

For black holes in GR, the analog of second law of ther- 
modynamics is the "area theorem" which asserts that 
area of a black hole can not decrease in any classical 
process A simple illustration of this result can be 
obtained by studying the physical process version of the 
first law which describes the dynamical change of horizon 
area in response to a flux of matter through the horizon 
Q. In case of stationary bifurcate horizons with a hori- 
zon generating Killing field £ a , the physical process first 
law is, 
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where k is the surface gravity of the black hole, T a b is the 
stress energy tensor of the matter and the integration is 
over the horizon H. The measure is dH b = —kb dX^j dA, 
where A is an affinc parameter along the horizon genera- 
tors, k b = (d\) b is tangent to the generators, and y^ydA 
is the area element of a constant A horizon slice. ( For a 
detail discussion and derivation of this law for GR, see 
Refs. 0, @ and @). The main ingredients in the deriva- 
tion are the Raychaudhuri equation and the assumption 
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of sufficiently quasi stationary process in which the ex- 
pansion and shear of the generators are small and hence 
all second order terms can be neglected. There is also 
an assumption that the black hole is stable under small 
perturbations. Eq.([T]) shows that in a dynamical process, 
the area of a black hole can not decrease as long as the 
matter satisfies null energy condition (i.e T a bk a k b > 0). 
Once we identify the Hawking temperature as T = k/2tt, 
Eq. (pj becomes completely analogous to the Clausius re- 
lationship AS = AQ/T as in ordinary thermodynamic 
systems. 

The question of validity of the full second law for ar- 
bitrary theory of gravity still remains an unresolved is- 
sue. Except for the case of /(i?)-gravity there is 
no proof of the analog of Hawking's area theorem be- 
yond GR. In the quasi-stationary case, an argument for 
second law valid for all diffeomorphism invariant gravity 
theories was given in Ref. [loj . under the assumption 
that the stationary comparison version of the first law 
implies the physical process version for quasi-stationary 
processes. Still, we need a direct proof of the physical 
process version, by computing the change of the Wald 
entropy due to the accretion of matter by the black hole. 
This is an important consistency check for the validity of 
black hole thermodynamics beyond GR [ll[ . If the phys- 
ical process and the equilibrium state versions of the first 
law do not match, this will suggest an inconsistency in 
the assumptions behind the derivation, in particular to 
the assumption that the black hole horizon is not de- 
stroyed by throwing matter into it and the final state in 
the late time is again stationary. The fact that in GR, 
both these versions agree, lends strong support to the 
idea that gravitational collapse result in a predictable 
black hole. It is expected that black holes in any reason- 
able theory of gravity must have such property. Also, a 
direct proof of physical process version of first law will 
automatically establish a quasi-stationary version of the 
second law [10l ). 

In this letter, we provide a direct proof of the physical 
process version of the first law for Einstein-Gauss-Bonnct 
(EGB) gravity. Our proof involves the direct evaluation 
of the change of the Wald entropy when a stationary 
black hole is weakly perturbed and finally settles down 
to a new stationary state. We show that in such a pro- 
cess, the Wald entropy for EGB theory can not decrease 
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as long as the matter stress energy tensor satisfies null 
energy condition. This firmly establishes that the Wald 
entropy for stationary black holes merits the name en- 
tropy even for more general gravity theories. 
The letter is organized as follows: we begin by presenting 
the properties of stationary Killing horizons and discuss 
how Einstein tensor takes a highly symmetric form on 
the Killing horizon of any stationary (and non-extremal) 
black hole spacetime Next, we review the Wald en- 
tropy for stationary horizons in EGB theory and present 
the proof of the physical process version of the first law. 
Finally, we conclude with some discussions 1 . 
Let us begin with properties of a stationary black hole 
horizon. In a D-dimensional spacetime, the event horizon 
is a null hyper-surface % parametrized by an affine pa- 
rameter A. The vector field k a = (d\) a is tangent to the 
horizon and obeys geodesic equation. Consider a spatial 
cross section at A = everywhere. All such A = con- 
stant slices are space-like and foliate the horizon. Any 
point p on such slices have coordinates {X,x A } where 
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, D) are the coordinates of a point on 



A = slice connected with p by a horizon generator. We 
can construct a basis with the vector fields, {k a ,l a ,e'X} 
where l a is a second null vector such that l a k a — —1. The 
induced metric on any slice is j a b = g a b + 2/c( a 4)- The 
change of the induced metric from one slice to another 
can be obtained from the metric evolution equation [l2j , 
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where a a b is the shear and 6 is the expansion of the hori- 
zon. If the event horizon is also a Killing horizon 2 , i.e. 
the horizon generators are the orbits of a Killing field 
£ a = (d/dv) a , which is null on the horizon, we can define 
the surface gravity k by the relationship £ a V a £ b = 
For stationary spacetimes with a Killing horizon, both 
the expansion and shear vanish and using Raychaudhuri 
equation and the evolution equation for shear, we obtain 
[H, R ab k a k b = fc a fc c 7™7nCa&cd = where C abcd is the 
Weyl tensor. 

The Einstein tensor (and hence the Ricci tensor) takes 
a highly symmetric, block diagonal form on a station- 
ary Killing horizon with a regular bifurcation surface 
[l3[. To understand this result, let us construct a basis 
{£ a , N a , e\}, where N a is another null vector satisfying 
£, a N a = — 1 and express Ricci tensor in this basis. On 
the stationary Killing horizon, Ricci tensor satisfies the 
relation R ab £, b ~ £a 12j. Using this, it is straightfor- 
ward to show, that on the bifurcation surface, the Ricci 



1 We adopt the metric signature (— , +,+,+,...) and our sign con- 
ventions are same as those of 1131 . 

2 Here we make an implicit assumption, that the event horizon 
of a stationary black hole is also a Killing horizon with regular 
bifurcation surface. Although this is certainly true for GR [l4j . 
we are not aware of any proof for EGB gravity. 



tensor must have a form, R a b — Cig^ b + C ab , where 
g^ = — 2fc( a Z b ), is the metric of the 2-dimensional plane 
orthogonal to the (D— 2)-dimcnsional bifurcation surface. 
The coefficient C\ is given by, C\ = —R a bk a l b and C ab is 
a tensor entirely intrinsic to the horizon cross section such 
that k a C ab = l a C ao = 0. Since Ricci tensor is a Killing 
invariant, we can Lie propagate the Ricci tensor onto any 
other space-like section of the horizon (For a related dis- 
For stationary Killing horizons, since 
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both the expansion and shear vanish, using the block di- 
agonal form of the Ricci tensor on the horizon, it is also 
possible to show that [HI, k a ^^ C jlk^a,bcd = 0. 
Another relationship which will be useful for our cal- 
culation is the expression of intrinsic Ricci curvature 
(D~2)j^ab f -j- ne cross _ sec t;i ori f the stationary Killing 

horizon in terms of full curvatures, given by [17| 



(D-2) 



Rab = la lb Is R-n 



(3) 



Next, we discuss the features of EGB gravity theory. A 
simple modification of the Einstcin-Hilbcrt action is to 
include the higher order curvature terms preserving the 
diffcomorphism invariance and still leading to an equa- 
tion of motion containing no more than second order time 
derivatives. In fact this generalization is unique (Tc| and 
the lowest order correction appears as the Gauss-Bonnet 
(GB) term, C GB = R 2 - 4:R ab R ab + R a bcdR abcd in space- 
time dimensions D > 4. EGB gravity is free from ghosts 
[ljl H(| , leads to a well-defined initial value problem and 
therefore, is a reasonable candidate for a low energy ef- 
fective theory of gravity. The action functional is given 
by, 
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The field equation of EGB theory is, G a b + ceH a b = 87rT b 
where, 
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For the action in Eq.Q, the Wald entropy associated 
with a stationary Killing horizon is [2l| , 



S = \J ' pVldA, 
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where the entropy density p = (l + 2a ( D ~ 2 ^R) and the 
integration is over (D — 2)-dimcnsional space-like cross- 
section of the horizon, ^"^i? is the intrinsic Ricci scalar 
of the horizon cross-section. 

Our goal is to prove that this entropy always increases 
when a black hole is perturbed by a weak matter stress 
energy tensor of order 0(e) provided the matter obeys 
null energy condition. Since the black hole is stationary 
in the asymptotic future, the vector field £ a is an exact 
Killing vector at late times and all the Lie derivatives of 
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dynamical fields w.r.t £ a vanish. The accretion process is 
assumed to be slow such that all changes of the dynamical 
fields are first order and we can neglect all viscous effects. 
More specifically, we assume that: 
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The bifurcation surface at A = is taken as the initial 
cross-section and the final cross-section is in the asymp- 
totic stationary regime. Then the change is entropy is 
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In deriving the second line, a total derivative is discarded 
since it vanishes both on the initial A = slice and on the 
final stationary state. Also, we have used the Raychaud- 
huri equation and neglected all terms except those first 
order in perturbation. Using field equation, this entropy 
change is expressed as, 

AS = J \{2nT ab -^n a ^ k a k b dX^jdA, (9) 

where K ab = H ab - 2^ D ~^RR ab + 2 V a V & P-^R. The 
first term in Eq.Q is linear in perturbation. Our aim 
is to prove that the first order part of the term H a bk a k b 
vanishes identically. In order to sec that, we note that the 
terms in lZ ab involve squares of curvatures. We are only 
interested in quantities first order in perturbation over 
a background stationary spacetime. The background 
spacetime in our context has no particular physical mean- 
ing except as a reference spacetime. Therefore, when we 
encounter a product of two quantities X and Y, to ex- 
tract the part linear in perturbation, we will always ex- 
press such a product as, 
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where X^ is the value of the quantity X evaluated 
on the stationary background and X^ is the per- 
turbed value of X linear in perturbation. Note that, on 
the stationary background, Raychaudhuri equation de- 
mands R l -f b k a k b = and since T^ ] k a k b = 0, we have 
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H K ab 'k a k° = 0. Also, to simplify the calculation, we use 
diffeomorphism freedom to make the null geodesic gen- 
erators of the event horizon of the perturbed black hole 
coincide with the null geodesic generators of the back- 
ground stationary black hole As a result, the pertur- 
bation in the location of the horizon vanishes. 
As an illustration of the perturbation scheme mentioned 
in the Eq. (fT0|) . we evaluate a term in H ab k a k b as, 
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where we have used the results of [13| to express the back- 
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ground Ricci tensor on the horizon as R y ab ' = Cig^ b +C ab . 



Also, note that l l ° R (P b] k a k b = R\[ b >k a k b . Implementing 
similar scheme for other terms, the first order part TZ\ of 
(H ab - 2 { - D -^RR ab ) k a k b becomes, 
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The last term in lZ ab k a k b requires special care. To eval- 
uate this term, we note that the change of the (D — 2)- 
dimensional scalar curvature can be thought of due to 
the change in the intrinsic metric. Then, we can calcu- 
late this change by using the standard result of variation 
of Ricci scalar as, 
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where D a is the covariant derivative intrinsic to the hori- 
zon cross-section and since the sections of the horizon are 
compact surfaces without boundaries, the surface term 
D a (S\V a ) in Eq. (fT3"| does not contribute. Then using 
Eq.([2]) and neglecting terms of higher order, we obtain, 
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where we have again used Raychaudhuri equation and the 
evolution equation for shear [l2[ keeping terms linear in 
perturbation. The last line follows after expressing the 
perturbed Weyl tensor in terms of curvature and Ricci 
tensors. 

We rewrite Eq. (fT4]) using Eq.Q and express R S mm in 
terms of Weyl and Ricci tensors of the background to 
arrive at, 
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where we have defined, C = 7 a& C a f,. Next, using the 
expression, R^ = 2Ci+C, and comparing with Eq. (fT2j) . 
it is straightforward to prove that the first order term in 
TZ ab k a k b vanishes identically and we have the result, 



AS = 
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+ 0(e 2 ). (16) 
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For the background stationary horizon, the Killing vector 
£ a is related with the horizon generators k a as £ a = nXk a , 
since T a \, itself is of 0(e), we can use this relation in 
Eq. (|16p . Then, identifying the Hawking temperature as 
T = k/27t, in the leading order of perturbation, we finally 
obtain, 

TAS=- [ T ab £ a dH b . (17) 

This is the desired form which establishes the physical 
process version of the first law for EGB gravity. If the 
matter stress tensor obeys null energy condition, Eq.(|16[) 
shows that the Wald entropy for stationary black holes 
in EGB theory can not decrease in a dynamical process 
which perturbs the black hole and leads to a new sta- 
tionary state. Since, it is expected that classical matter 
obeys null energy condition, we can conclude that as in 
case of GR, the entropy of stationary black holes in EGB 
theory can not decrease through any classical process. 
An interesting feature of this derivation of physical pro- 
cess first law is that it is completely local in nature. Un- 
like the equilibrium state version, there is no reference of 
the asymptotic infinity. In fact, it is the field equation 
which enforces the first law for quasi stationary changes 
of the horizon. In case of GR, it is possible to reverse this 
argument and derive Einstein equation as an equation of 
state of the spacetime [22| ■ Attempting similar construc- 
tion for EGB gravity will require a notion of entropy for 
non-stationary states. For some possible proposals, see 
Refs. [i| and |j. 

The derivation of the laws of black hole mechanics is en- 
tirely classical. To complete the thermodynamic analogy, 
we need to invoke quantum theory. The idea that black 
holes radiate at a temperature T = k/2it is entirely a 
consequence of quantum field theory in the presence of 
a horizon and independent of the gravitational dynam- 
ics. As a result, the mere analogy of classical black hole 
physics with thermodynamics becomes an exact law once 
the quantum effects are taken into consideration. Then, 
a natural interpretation of the black hole entropy is that 
it counts the quantum micro-states of the black hole. For 
any reasonable theory of gravity, which has stable black 
holes, the density of states must be of the form exp (Sw), 
where Sw is the Wald entropy. In fact, in the context 
of string theory, it has been shown at least for extremal 
and near-extremal black holes in EGB theory, that the 
microscopic computations exactly matches with the Wald 
formula [2^,|2{|. Also, the GB term appears as a low en- 
ergy a' correction in case of tree level heterotic string 
theory (l9| . Hence, from this microscopic point of view, 
it is desirable that the Wald entropy of the black holes 



in EGB theory should increase and our result establishes 
this at least in the context of first order perturbation. 
For ordinary thermodynamic systems, the entropy, by 
construction is a path independent function and the 
change of the entropy from one stationary state to an- 
other is always given by the difference of entropies be- 
tween two states, independent of the process. The agree- 
ment of the equilibrium state and the physical process 
version of first law establishes the same property for black 
hole entropy. Hence, our result implies that the thermo- 
dynamic analogy for black holes is not only a feature of 
GR and simple f(R) theories, but equally applicable for 
well-motivated modifications of GR, like the EGB the- 
ory. For f(R) gravity, it is not very surprising that the 
physical process version holds, since all f(R) theories are 
related to GR by a field dependent conformal rescaling of 
the metric. But, the validity of physical process version 
of first law for EGB gravity is a non-trivial check suggest- 
ing a more general nature of black hole thermodynamics. 
A possible generalization of our work would be to study 
black holes in a general Lovelock theory [l8[ . For any ra- 
th order Lovelock term C m in D-spacctime dimensions, 
the entropy density of the stationary black holes is (2~i| 
simply proportional to ( D ~ 2 '£r m -i), the previous Love- 
lock term but intrinsic to the horizon cross-section. This 
remarkable property suggests that our method can be 
readily used to investigate the physical process version of 
the first law for any Lovelock theory. 
In this letter, we have only considered the increase of 
the black hole entropy alone, neglecting the contribution 
of the quantum fields outside the horizon. In GR, we 
know that the effective negative energy flux from quan- 
tum fields to a black hole can lead to a decrease of horizon 
entropy. Then the relevant question is whether a general- 
ized second law (A(Sbh + Soutsidc) > 0) holds. There are 
arguments (2?| that the generalized Second Law applies 
for semi-classical processes in case of black holes in GR. 
Validity of these arguments for reasonable higher curva- 
ture gravity theories is still an open problem. 
Finally, it is worthwhile to study whether the physical 
process first law holds for any diffeomorphism invariant 
gravity theory or applies to a special class of action func- 
tional. The answer, in either way, will be an important 
input for any quantum theory of gravity. 
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